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We study the massless scalar field theory in a four-dimensional Euclidean space, where all 
but one of the coordinates are unbonnded. We are considering Dirichlet boundary conditions in 
two hyperplanes, breaking the translation invariance of the system. We show how to implement 
the perturbative renormalization up to two-loop level of the theory. First, analyzing the full two 
and four-point functions at the one-loop level, we shown that the bulk counterterms are sufficient 
to render the theory finite. Meanwhile, at the two-loop level, we have to introduce also surface 
counterterms in the bare lagrangian in order to make finite the full two and also four-point Schwinger 
functions. 
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I. INTRODUCTION 

In this paper we are interested to show how to implement the renormalization procedure up to two-loop level in the 
massless scalar field theory, defined in a four-dimensional Euclidean space with one compactified dimension. Our 
aim is to shed light on the renormalization procedure in a system defined in a domain where translational symmetry 
is broken, which must be done for example in the high temperature dimensional reduced quantum chromodynamics 
(QCD), defined in a finite region. 

Quantum chromodynamics is a non-abelian Yang-Mills theory with gauge group SU(S). Since it is assumed that 
the fermions of the theory transform according to the fundamental representation of the gauge group, each flavour 
of quark is a triplet of the color group SU{3). Gauge bosons transform according to the adjoint representation. The 
interaction between the quarks is mediated by the gluons. Due to the non-abelian structure of the theory, the gluons 
conple not only with the quarks but have also cubic and quartic self-interaction. The self-interaction of the gluons 
provides the anti-screening of the color charge in QCD. This is responsible for asymptotic freedom and presumably 
confinement. 

The confinement-deconfinement phase transition in QCD may occur in usual matter at sufficiently high temperature 
or if is strongly compressed QQIl In ultra-relativistic heavy ion collisions, we expect that the plasma of quarks and 
gluons can be produced. We would like to stress that, although non-equilibrium processes occur in the quark-gluon 
plasma in the heavy ion collisions, for simplicity in a first approximation we can assume a static situation. Just after 
the collision hot and compressed nuclear matter is confined in a small region of the space and in such circumstances 
the volume and surface effects become very important. 

In the above described physical situation there are two important points: the first one is that the thermodynamic 
limit of the infinite volume system can not be used and therefore finite volume effects should be investigated and taken 
into account. The second point is that the quark-gluon plasma exists in a situation of high temperature, where using 
the Matsubara formalism to describe high temperature QGD, dimensional reduction must occurs i i i Q Dimen¬ 
sional reduction is based in the Appelquist-Garrazone decoupling theorem @ . From a more fundamental theory, the 
effective Lagrangian density of this theory can be obtained as some low-energy limit of the fundamental theory where 
the heavy modes have been removed. There are some interesting physical situations where the decoupling theorem 
can be used. First, for scalar fields without spontaneous symmetry breaking. Second, in quantum electrodynamics, 
where a derivative expansion of the photon effective action can be obtained by integrating out the fermionic fields. 
Also in QCD at least in lowest order in perturbation theory the decoupling theorem works. The decoupling theorem 
is not valid for example in spontaneous broken gauge theories. It is important to stress that the non-validity of the 
decoupling theorem means that low-energy experiments can provide information about the high energy physics. 

Going back to the heavy-ions collision situation, we can assume that the following scenario appears: in the situation 
where dimensional reduction occurs, we have an effective theory for the gluons field and also finite size effects for 
these bosonic fields. To shed light on the renormalization procedure in systems defined in domain where translational 
symmetry is broken, as for example the high temperature dimensional reduced QCD, in this paper we are interested 
to investigate scalar models, impose classical boundary condition over the fields. We hope that this study will give us 
some insight over the most interesting and also more complicated situation as the one mentioned above. Therefore, 
in this paper we analyze how to implement the perturbative renormalization up to two-loop level of massless 
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scalar field model defined in a four-dimensional Euclidean space with one compactified dimension. 

Finite size effects and the presence of macroscopic structures in different field theory models has been extensively 
studied in the literature. The critical behavior of the 0{N) model in the presence of a surface was a target of intense 
investigations 0- The same 0{N) model was studied in two different geometries: the periodic cube and the cylinder 
along one dimension (the time) and finite and periodic in the (d— 1) remaining dimensions by Brezin and Zinn-Justin 
M- Finite size effects in QCD im and also in different field theory models also has been extensively studied in 
the literature. Assuming periodic or antiperiodic boundary conditions for bosonic and fermionic models respectively, 
the translation symmetry is maintained, and surface effects are avoided. Therefore, to avoid surface effects, quantum 
fields defined in manifolds with periodic or anti-periodic boundary conditions in the spatial section was preferred by 
many authors [Q. Nevertheless, the case of boundaries conditions that break the translational symmetry deserves 
our attention. 

In the case of hard boundary conditions as for example Dirichlet-Dirichlet(DD) or Neumann-Neumann (NN), the 
translational invariance is lacking. This fact makes the Feynman diagrams harder to compute that in unbounded 
space. Moreover the renormalization program is implemented in different way from unbounded or translational 
invariance systems since some surface divergence appear [T^ . For translational invariant systems, one can use the 
momentum space representation, which is a more convenient framework to analize the ultraviolet divergences of a 
theory. Translational invariance is preserved for momentum conservation conditions. For non-translational invariant 
systems a more convenient representation for the n-point Schwinger functions is a mixed momentum coordinate space. 

Fosco and Svaiter considered the anisotropic scalar model in a d-dimensional Euclidean space, where all but one of 
the coordinates are unbounded. Translational invariance along the bounded coordinate which lies in the interval [0, L] 
is broken because the choice of boundary condition chosen for the hyperplanes at z = 0 and z = L. Two different 
possibilities of boundary conditions was considered: (DD) and also (NN), and the renormalization of the two-point 
function was achieved in the one-loop approximation |l4| . Further the renormalization of the four-point function was 
achieved in the one-loop approximation by Caicedo and Svaiter M Finally Svaiter M studied the renormalization 
of the massless scalar field model in the one-loop approximation in finite size systems assuming that the system 
is in thermal equilibrium with a reservoir. Also, still studying surface, edge and corners effects, Rodrigues and Svaiter 
Ea analysed first the renormalized vacuum fluctuations associated with a massless real scalar field, confined in the 
interior of a rectangular infinitely long waveguide. A closed form of the analytic continuation of the local zeta function 
in the interior of the waveguide was obtained and a detailed study of the surface and edge divergences was presented. 
Next, these authors M studied the renormalized stress tensor associated with an electromagnetic field in the interior 
of a rectangular infinitely long waveguide. 

In this paper we will consider an interacting massless scalar model, in a four-dimensional Euclidean space, where 
the first three coordinates are unbounded and the last one lies in the interval [0, L], We analyze only DD boundary 
conditions. First we present an algebraic expression in coordinate space for the free propagator which let us to identify 
the divergences of the n-point Schwinger functions for the interacting theory. This algebraic expression agrees with 
the result obtained by Lukosz M We would like to stress that instead of assume hard boundary conditions, some 
authors assumed soft boundary conditions and also treated the boundary as a quantum mechanical object [^. Here, 
we prefer to keep hard classical boundary conditions. 

The organization of the paper is as follows: In section II we discuss the slab configurations, obtaining some important 
expressions for the free propagator in order to understand some procedures in the divergences identification. In section 
III the regularization program is implemented in the one-loop approximation. In section IV the regularization program 
is implemented in the two-loop approximation. Section V contains our conclusions. In appendix A, expression for the 
free propagator is introduced. Throughout this paper we use h = c = I. 


II. CLASSICAL BOUNDARY CONDITIONS AND SOME PROPERTIES OF THE FREE PROPAGATOR 


Let us consider a neutral scalar field with a {X(p'*^) self-interaction, defined in a d-dimensional Minkowski spacetime. 
The vacuum persistence functional is the generating functional of all vacuum expectation value of time-ordered 
products of the theory. The Euclidean field theory can be obtained by analytic continuation to imaginary time 
allowed by the positive energy condition for the relativistic field theory. In the Euclidean field theory, we have the 
Euclidean counterpart for the vacuum persistence functional, that is, the generating functional of complete Schwinger 
functions. The {X<p^)d Euclidean theory is defined by these Euclidean Green’s functions. The Euclidean generating 
functional Z(h) is formally defined by the following functional integral: 


Z{h) = J [dip] exp S'o — Sj + J df‘xh{x)ip{x^ 


(1) 
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where the action that describes a free scalar field is 

5'oM = J d'^x , (2) 

and the interacting part, defined by the non-Gaussian contribution, is 

Sii^) = J d‘^x^ ip*{x). (3) 

In Eq.lpJ, [dt/j] is a translational invariant measure, formally given by [dip] — Y\^dp{x). The terms A and mg are 
respectivelly the bare coupling constant and mass squared of the model. Finally, h{x) is a smooth function that we 
introduce to generate the Schwinger functions of the theory by means of functional derivatives. Note that we are 
using the same notation for functionals and functions, for example Z{h) instead the usual notation Z[h]. 

In the weak-coupling perturbative expansion, we perform a formal perturbative expansion with respect to the non- 
Gaussian terms of the action. As a consequence of this formal expansion, all the n-point unrenormalized Schwinger 
functions are expressed in a powers series of the bare coupling constant qq. Let us summarize how to perform the 
weak-coupling perturbative expansion in the {\p'^)d theory. The Gaussian functional integral ZQ{h) associated with 
the Euclidean generating functional Z{h) is 

Zo{h) = -^ J [d‘f] exp ^-^pKp + hp^ . (4) 

We are using a compact notation and the first term in the right-hand side of Eq. 0 is given by 

pKp = j (fx J d!^yp{x)K{mQ\x,y)p{y). (5) 


The term that couples linearly the field with the external source is 

hp = J d^x p{x)h{x). ( 6 ) 

As usual A/" is a normalization factor and the symmetric kernel K{mp,x,y) is defined by 

Ar(TOo;a:,2/) = (-A -b mg ) - J/), (7) 

where A denotes the Laplacian in the Euclidean space As usual, the normalization factor is defined using the 

condition Zg(/i)|?i=g = 1. Therefore N = (det(—A -b rrig)) ^ but, in the following, we are absorbing this normalization 
factor in the functional measure. It is convenient to introduce the inverse kernel, i.e. the free two-point Schwinger 
function Go{mo;x — y) which satisfies the identity 

[ d'^z Gg(mg; x - z)K{mo; z - y) = 5'^{x - y). (8) 


Since Eq. 0 is a Gaussian functional integral, simple manipulations, performing only Gaussian integrals, gives 


J [dp] 



d‘^x 


d‘^y h{x) (mg; a; - y)h{y) 


(9) 


(‘ 2 '\ 

Therefore, we have an expression for Zo{h) in terms of the inverse kernel Gg ' (mg; x — y), i.e., in terms of the free two- 
point Schwinger function. This construction is fundamental to perform the weak-coupling perturbative expansion with 
the Feynman diagramatic representation of the perturbative series. The non-Gaussian contribution is a perturbation 
with regard to the remaining terms of the action. It is important to point out that the weak-coupling perturbative 
expansion can be defined in arbitrary gemetries, and classical boundary conditions must be implemented in the 
two-point Schwinger function. Another way is to restrict the space of functions that appear in the functional integral. 

We are interested to study finite size systems, where the translational invariance is broken. In this situation, we are 
analyzing the perturbative renormalization for the ^p‘^ massless scalar field model, in the two-loop approximation. 
Therefore, let us assume boundary conditions over the plates for the massless field p{x). For simplicity we are 
assuming Dirichlet-Dirichlet boundary conditions i.e.. 


p{r, z)[^=Q = p{r, z)[^=L = 0 , 


(10) 
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for the free field. Since the translational invariance is not preserved, let us use a Fourier expansion of the fields in the 
following form 


(p{f, z) = 


( 2 ^)- 




( 11 ) 


where the set u„(z) are the orthonormalized eigenfunctions associated to the operator — 

,2 

and kn = n = 1,2.... The orthonormal set correspond to the eigenfunctions of the Hermitian operator — ^ defined 
on a finite interval is given by 


, , /2 . ,mrz 

Un{z) = Y ^ n = l, 2 ... 


These eigenfunctions satisfy the completeness and orthonormality relations, i.e.. 


' Un{z)ul{z') =S{Z- Z') 


( 12 ) 


(13) 


and 


dzUn{z)u*^,{z) = 6n,r 


(14) 


Since we are interested to perform the weak coupling expansion, let us first write the free two-point Schwinger function. 
This free two-point Schwinger function can be expressed in the following form 


G''o\r,z,z') = 


1 


(27r) 


d-1 


d‘^ ^p^e^P-^Un{z)u*^{z')Go,nil^, 


(15) 


where Go^nip) is given by 


Go,n{P) = (P^ + + 


(16) 


Next, we will present some properties of the two-point free Schwinger-function in order to understand the behavior 
of the interacting field theory in the presence of macroscopic structures. Therefore, in order to understand some 
procedures used in the identification of the divergences in the Schwinger functions that will appear in the next 
section, let us analyze some properties of the freet two-point Schwinger function. Substituting Eq. m and Eq. fTT)|) 
in Eq. m we get that the free propagator G^\fi — r 2 , zi, Z 2 ) can be written as 


00 


(ri - r 2 , zi, Z 2 ) = ^ 2^ sm 


. ^mrz2^ 


dd-lp gip.(Fl-f2) 

(27r)^“l _|_ (1^)2 _|_ ^2^ 


(17) 


The next step is to show that the two-point free Schwinger function can be written in terms of the variables: ri 2 , Z 12 
and finally 2 ^^, where ri 2 = ^12 — respectively. Working in the four-dimensional case 

and also in the massless situation, a straightforward calculation (see appendix A) give us that G^'^ (rq — fq, Zi, Z 2 ) can 
be written as: 


G'o^^(ri - r 2 ,zi,Z 2 ) = 


1 




E 

k— — oc) 


(fc-^) 


^2 + (1^)2 


{k- 


L)2 + (112)2 


(18) 


The former expression for the two-point Schwinger function was obtained also by Lukosz |l9| using the image method. 
Performing the summations in Eq. m (see appendix A), it is possible to find a closed expression for GQ^^(ri—fq, zi, Z 2 ). 
We get 


G''o\ri - r2,zi,Z2) 


sinh(7rri2) 

sin(2£i)sin(^) 

IGttL^ ri2 


Sinh2(in2)+sin2(^) 

sinh2(^) -hsin2(ll^) 

> 


(19) 
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It is not difRcult to show that the two-point Schwinger function Gq '^(ri — ^ 2 , 21 , 22 ) satisfies the following properties: 

(o'] 

(i) The free two-point Schwinger function is not negative, i.e., Gg '{fi — 7 ^ 2 , 21 ,^ 2 ) > 0 , for zi, Z 2 G [0, L] and 
ri,r2 G since we are working in an Euclidean space. 

(ii) The free two-point Schwinger function is zero when one of its points are evaluated on the boundaries 

in -^2,0, Z2) = G^^^ (fi - r2, L, Z2) = (n - r2,zi, 0 ) = g[,^^ (ri - r2, Zi, L) = 0 , 


since we are assuming Dirichlet boundary conditions. 

(iii) The free two-point Schwinger function contain the usual bulk divergences, i.e., when (ri,zi) = (T 2 ,Z 2 ), it is 
singular. From the Eq. ra we can identify three singular terms. Splitting the free two-point Schwinger function in 
the singular and regular terms we have: 


Go\fi-f2,Zi,Z2) = 


1 


1 


47r2L2 I )2 

1 I 1 

'47r2L2 


(2 - zJ * 2)2 + r ^2 


00 

E 


rl, (2fc- 


-12I 


12 


~ ^ 12 )^ 


( 20 ) 




k 

k^0,l 


12 


The first term of the right side of the last equation, is singular only when rq = fq and zi = Z 2 - This is the term that 
carries the usual bulk divergences. The second term is singular only when zi = Z 2 = 0 and rq = rq. The third term 
is singular only when z\ = Z 2 = and iq = fq. These two terms mentioned previously, carries surface divergences. 
Finally the two last terms do not have singularities. 

{iv) When |rq — rql/L » 1 the free propagator behaves like 


G^i^iri -r2,zi,Z2) 


27rL2 


ri2 


■ Sin 


(^) 


Sin 


iT 




( 21 ) 


which shows an exponential convergence behavior. 

(v) The integral of the variable {r,z} on a neighborhood around {r',z'} of the free propagator is finite, i.e., 
(Prdz G^Q^f— r', z, z') < 00 . See Fig.l^). 

Property (v) allow us to shows that the external legs of the Feynman diagrams do not create divergences. Let us 
suppose we have the integral corresponding to some Feynman diagram, 


[ (frdz G^o\r-r',z,z') F{r',z') 
Jr 


( 22 ) 


where Gg {r — r', z, z') is some external leg and F{r', z") describes the remainder part of the diagram. Now in order 
to proceed we have to use the following statement: for two continuous and positives functions /(a:) and gix) defined 
in a finite region i? with the exception of the point x\ where /(af) diverges, then the integral / = f{x)g{x) is 

finite, if and only if F = jy d‘^x f{x) is finite on some neighborhoods V of the point xi. With the property (v) and 
the statement before we can see that external legs from the Feynman diagrams do not generate divergences. 


(f,2) 



FIG. 1: 


III. REGULARIZED TWO AND FOUR-POINT SCHWINGER FUNCTIONS AT ONE-LOOP ORDER 

In this section we identify the divergent contribution in the two and four-point Schwinger function at one-loop level. 
Essentially we use the Eq. in the IPI diagrams of the Green functions considering their external legs, and the 
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integrations in the coordinate space. We write Eq. Pl|l 


as 


G''o\fi - r2,zi,Z2) = 


where h{ri 2 , zi, Z 2 ) is given by 


l(zi2)^+r?2 (^i)^+^i2 (2-z^2)^+r^2 


h{ri2,Zi,Z2) 


(23) 


h{ri2,zi,Z2) = - 


4 .frl (fc- ^)2 + (m)2 




fc^O, 1 


(fc- fia)2 + (m)2 


(24) 


From the property {Hi) we see that the three first contributions in the right hand side of Eq. PI have singularities. 
Otherwise, the last term is finite in the whole domain where we defined the model. After this briefly introduction, we 



FIG. 2: The two-point function at one-loop level 

are able to study the interacting theory. Let us start analyzing the tadpole diagram, displayed in Fig. from which 
we can write the expression for the one loop two-point Schwinger function G) {ri — f 2 , zi, Z 2 ). We have that 


Gf ^ (ri - r 2 , zi, 22 ) = ^ J d^rdz g[,^^ (fq - f, zi,z) G^,^^ (0, z, z) G^,^^ (fq - f, Z 2 , z) . 


(25) 


In the following we are generalizing the results obtained by Fosco and Svaiter M- Let us begin studying the quantity 
Go^^(0,z,z) that appear in the tadpole defined in Ea. pil . From Ea. ll^ we get that Go^^(0,z,z) can be written as 


G'i\o,z,z) = 


47r2L2 


A- 


1 


1 


00 CXD 


1 


(2z/L)2 (2 - 2z/L)2 (2fc)2 (2k - 2z/L)2 

fe^O, 1 


k = - 
feT^O 


(26) 


where A is given by 


A = 


lim 


L2 


(zi.ri)^(z2T2) (Zi - Z2)2 -I- |fi - fqp 
A2 , 


L^S4 


= lim 

A->oo 87r2 


(27) 


Sd = r(d/2) ultraviolet cutoff. In the same way, from En . (12611 by performing the summations, we get for 

G^o"\0,z,z) 


g(')(0,z,z) = 


Substituting Eq. (HZI) in Eq. we obtain 


1 


47r2L2 


A 

- 


1 


12 4 sim(7rz/L) 


^4 


G^o^\0,z,z) = hm Tr^A^ 


A^oo 327r4 48L2 167.2 sin2(7rz/L)' 


(28) 


(29) 
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The first term in Eq. is a bulk divergence. Substituting Eq. 12911 in Eg. 12511 we get 


G^i\n -r2,zi,Z2) 


+ -r,Z 2 , z) 

^ /" j3 j (^1 -r,zi,z) if2 -f,Z2,z) 

-• 


(30) 


The first term in the right hand side carries a bulk divergence. The second term is finite. To see this we analyze the 
integral by sectors. Therefore we have 



f dPrdz (fi -f,zi,z) G'^q^^ ir2 - f, Z2, z) = f 
Jr Jri 



(31) 


where each integral are defined in different regions displayed in FigEl where the points (ri,zi) and (^ 2 ,^ 2 ) are the 
centers of the regions Ri and R 2 respectively. Using the property (v) we have that the integrals on i?i and R 2 
are finite. Since the free propagators G'Q^^(ri — f,zi,z) and G'^\f 2 — f,Z 2 ,z) presented in Eq. ED does not have 
divergences on R 3 and this region is compact, then the integral on R 3 is finite. The integrals defined in regions i ?4 
and i ?5 also are finite since from the property (iv) the propagator decreases exponentially when one of its points 
become far from the other. Thus the integral defined by Eq. m is finite. Finally we have to study the third integral 
in the right hand side of Eq. CT . Note that the term diverges when z is evaluated on the boundaries. 

/n'i ^2^ 

Nevertheless this integral is convergent, because the products of Gq '(ri — r, zi, z) and Gq '(r-j — r, Z 2 , z) take away 
the divergence. Using Eq. m, we have that third integral in the right hand side of Eq. is finite. Therefore the 
one-loop two point Schwinger function only has bulk divergence. 

Our next step is to analyze the four-point Schwinger function in the one-loop level. Since the free propagator only 
has singularities when its two points are equal or also when this two points joined are evaluated at the boundaries, 
we continue our analysis of the integrals only in the domains where the two external points of the free propagators 
take the same values. The complete four-point function at one-loop level is given by 


G^‘‘^(fi,zi,f2,Z2,r3,Z3,f4,Z4) = Y y J j rf-z'- f, zi,z)G[,^\f 2 - r, Z 2 ,z) 


G^g\r- r',z,z') G{,'^^(r 3 - r',Z3,z') - r",Z4,z'). 


.(2) 


( 2 ), 


(32) 


For simplicity, in Fig. we define three different regions between the boundaries. The first one, i?i is concerned when 

{r', z'} is close to {r, z}. In this region the contribution coming from Gg ’ {r — r', z, z') is singular. Nevertheless, 

we still have to analyze if this divergent behavior will appear in the integral defined by Eq. ED- We will show 
that the singularities will appear only as bulk divergences. In the region R 2 (z, z' —> 0 and r'^r) the term 


G[,")(r-^ 


') is also divergent. As we will see, this divergent behavior disappears when we compute the 


















FIG. 4: The four point function at one loop. 



FIG. 5: Regions of integration for the four point function 


complete four-point function at one-loop order, defined by Eq. 123). In the region i ?3 (z, z' ^ L and r' ^ the 
situation is identical as in the region i? 2 . Using the same argument that we used before to analyze the convergence of 
the integral defined by Eq. 123, we can study the convergence of the integral defined by Eq. 123 with the amputated 
external legs. Therefore we have to study Eq. 123 with the external legs amputated. Therefore we have to study the 

r ( 2 ) 1 ^ 

quantity J d^rdzd^r'dz' Gq {f'— r, z', z) . Substituting Eq.(|23 in the former equation we get: 


d^rdzd^r'dz' 


-r,z',z) 


1 


(dTT^L^)' 


■ (^h + h + h + h + + finite part 


where the integrals li, i = 1,2,.. are given by 

1 


h = 

h = 

h = 
h = 
h = 


d^rdzd^r'dz' 

d^rdzd^r'dz' 

d^rdzd^r'dz' 

d^rdzd^r'dz' 

d^rdzd^r'dz' 


[(^ 12 )^ + ^ 12 ] 

1 


[(z+)2 -h rf^] [(2 - z+)2 -h r22]' 


[(2=12)^ + ^? 2 ] [(^2)^ + ^^2] [(2=12)^ + rh ] [(2 - 42 ^ + 4 ] 

2 


[( 2 +)^ + 4] [(2- 42V + 4] ’ 

1 1 


1 


.( U 2)^+^12 ( 2 :^ 2 )^ +^12 ( 2 -^:+ 2 )^+^ 12 - 


h{ri2,zi,Z2 ). 


(33) 

(34) 

(35) 

(36) 

(37) 

(38) 


Let us investigate each term of Eq. 123- The integral /i must be analyzed only in the region Ri. For this purpose we 
need an auxiliary result. We can prove that a continuous and positive function f{x) which does not have singularities 

except for a: = 0, and M = where w'^ = |w;p, then there exist e' such that M = Sd Jq dw f {w"^) 

where e < e' < Vd e. Then we get 


rr+e 


It = 


d^r'dz' 




iRi 


[42^+4 


2 12 


z+e 


dz' 


[(z-U)2 + |r-r'|2]2 
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dw —j = S '4 In w 


(39) 


Therefore Ii contributes with a bulk divergence of the type as the one that appears in the theory without boundaries. 
In the usual renormalization procedure, the contribution coming from Ji can be eliminated by the usual counterterms. 
Concerning the contribution coming from I 2 we have that the first term +^ 12 ]^ singular in the region 

i?i. In the region i? 2 , using the same auxiliary result that we used before, we can obtain an upper bound to the 
contribution coming from this term. We get 



rdzd^r'dz' 


d^r 


[{z + z'Y + \r 

d w — = —b'ae 
:• 12 




pr+e /.£ /.£ 

dV / d^r' 


dzdz' 


0 Jo 


+ r - r' 


(40) 


Since the region R' C R 2 is finite this integral is convergent. Next, let us analyze the term l/[(2 — z(^)^ + ri 2 ]^ of I 2 
in the region R 3 . Since the behavior of the field in each plates (for z = 0 and z = L) are the same, then the analysis 
follows the same lines as previous ones and therefore this contribution is also finite. To study I 3 , we consider first the 
term 2 /[(zC 2 )^ + rj 2 ][(zi 2 )^ + ^ 12 ]- This expression must be studied in the regions i?i and i ?2 respectively. In Ri we 
can see that the convergence of 


f d^rdzd^r'dz' --7-5- , 

/fli [(z - z'Y + |r - r'p] [{z + z )^ + |r - f p] 

'-V-' 

finite in Ri 


(41) 


depend of the convergence of 


d^rdzd^r'dz' 


1 


{z- z')'^ + \r-f'\ 
From above arguments we have that Eq. 1221 ) can be written as 


d'^w 




= 84 


dw w = 


S' 4 e' 


/ 2 


(42) 


(43) 


thus Eg. 14211 gives a finite contribution. Now we consider the first term of I 3 in the region i? 2 . For this purpose we 
will use the following property. Let us take a continuous and positive function f{x) which does not have singularities 

except for a; = 0, and N = /J d^yd^z f{y‘^ + z^) then there exist e' in such way N = dw /(w^) 

where d > 0. Using this property, we have for the first term of Eq. (ISSll . in the region i? 2 , that 


nz-\-e 


d^r I dz I dz' I d^r' 


fo 


< / d^r / dz du d^v 
J Jo Jo J-e 


[(z - z')2 + |F_ r'P][(z + z')2 + |F- r'P] 


< -S '4 / d'^r J dz j dw ■ 


(u^ + v^)(z^ + + v^) 4 y Jo Jo (z2 + u;2) 25 'i52 


t". (44) 


Therefore the first term of lo is also finite in i? 2 - The second term 2/[((z(^)^ + — ^ 12 )^ + ^ 12 )] -^3 must be 

analyzed also in the regions Ri and R 3 . This analysis follows the same lines as the last case, therefore the contribution 
coming from this term is also finite. 

We have now to study the term I4. Note that 2/[((zj^)^ + — ^ 12 )^ + ^ 12 )] must be analyzed in the regions 

i ?2 and i ?3 respectively. Let us start with the region i? 2 - Using previous arguments we have that the convergence of 
I4 depends of the convergence of the following expression 


dzdz' 


^r+e 


/ 

a r - 


Si 


0 Jo 


z^ + z'^ + \r - r' 


= — dw — 7 - = — e 

4 Jo w^ 12 


5 / 3 


(45) 


which is finite. In the region Ro our analysis follow the same lines as in the region R 2 , thus the integral in the region 
i ?3 is also finite. 

Using the same argument that we used before, is not difficult to show that the contribution coming from I 5 is also 
finite. We conclude that the integral given by Eq. 021 ) only have bulk divergences. In this way we can conclude that 
at one-loop level the bulk counterterms are sufficient to render the complete connected Schwinger functions finite. In 
the next section we will identify the divergent contribution in the connected two-point Schwinger functions at the 
two-loop order. 
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FIG. 6: Two point Schwinger functions at two-loop level 


IV. THE DIVERGENCES IN THE TWO-POINT SCHWINGER FUNCTIONS AT TWO-LOOP LEVEL 


In this section we will generalize some results obtained by Fosco and Svaiter [Q and also by Caicedo and Svaiter 
M We will identify the divergent contribution in the connected two-point Schwinger functions at the two-loop order. 
The diagrams that we are interested to analyze are displayed in Fig. ®. The expression that corresponds to Fig. lOa) 
is given by: 


\ 2 /* r n 2 

^ ,^(2) l\ ^(2),^! ^ / X ^(2) 


^ / (fr'dz'd^rdzG''Q’{fi-f',zi,z') G}f\f'-f,z',z) {0, z, z) G';^\r 2 - r', Z 2 , z'). 


(46) 


Since the external legs in Eq. (gnj does not contribute to generate divergences, let us consider only the following 
integral 


d^rdz 


Replacing Eq. in Eq. lITTI) we get 
/ d^rdz 


G^o\r' -r,z',z)\ G^q\o,z,z). 


G‘'o\f'-f,z',z)\ G^q\0,z,z) = \im I d^rdz (r' - r, zz) 


(47) 


A^oo . 
1 


48L2 

1 

'l 6 L 2 


d^rdz 

d^rdz 


G^Q\r' - r,z',z) 


Go^\r' - r,z',z) 


sim(7rz/L) 


(48) 


The first term and the second one in Eq. gHJ can be renormalized introducing only bulk counterterms. The most 
interesting behavior appears in the last term of this equation. Note that this unrenormalized quantity contains only 


bulk divercences, since the contribution coming from 


G(^)(r'-f,z',z) 


1 2 


cancels the surface divergent behavior 


generated by the ^ term. Nevertheless, after the introduction of a bulk counterterm to render the contribution 

Gq {r' — r,z',z) finite between the plates, surface divergences appear. Thus this surface divergences must be 
renormalized. After the introduction of surface and bulk counterterms, the finite contribution coming from the last 
term of Eq. gHI), up to a finite renormalization constant, is given by 


1 


16L2 


d^rdz 


(Gf\r' 


-r,z,z) - 


1 


1 


1 


(47r2L2)2 _(_ ^2^J2J [sin^(7rz/L) (Trz)^ tt‘^{L-z)^_ 

Therefore this term contains an overlapping between bulk and surface couterterms. 

We still have to analyze the sunset diagram. The expression corresponding to the Fig. ®b) is given by 

\ 2 r 3 

— / d^r'dz'd^rdzG‘'Q\ri-f',Zi,z') (r' - r, z', z)J G^^^ (fa - r', Za, z')■ 


(49) 


(50) 


Again, the external legs does not contribute to generate divergences, and therefore let us study the amputated diagram, 
i.e., without external legs. We have 


d?rdzd^r'dz' 


Go\r' -r,z',z) 


(47r2L2) 


>2,2)3 (-^1 


+ I 2 + ... + /12 + finite part, 


(51) 
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where 


/i = y SrdzSr'dz' ■ 

I2 = J d^rdzd^rdz’ ■ 

h=jSri.4h'i-Z 

/4 = y d^rdzd^r'dz' ■ 

,,=Jd^rdM' 
/, = /dW.-V4 

/y = y d^rdzd^r'dz' ■ 
/s = y d^rdzd^r'dz' ■ 
/g = y d^rdzd^r'dz' ■ 
/lo = y d^rdzd^r' dz' 

111 = y d^rdzd^r'dz' 

1 12 = [ d^rdzd^r'dz' 


[{^i 2 ?+rl 2 ?' 

1 

1 

[(2-4)^ + ^? 2 P’ 

2 

2 

4r2)" + 4]^[(2-4)" + 4]’ 

2 

4r2)^ + 4][(4)^ + 4P’ 

2 

[fe)^ + 4] [(2 - 4)^ + 4]^ ’ 

2 

[(4)" + 4]^[(2-42)^ + 4]’ 

2 

[(4)^ + 4][(2-4)^ + 4F’ 

6 

[fc)^ + 4] [(4)^ + 4] [(2 - 4)^ + 4]' 

1 1 1 


-42)^+ 4 (4)^ + 4 (2-4)^ + 4- 

1 1 1 


(zi2)2 + 4 (4)^ + 4 (2-4)^+ 4 


h{'>'l2,Zi,Z2), 

h'^iri2,Zi,Z2). 


(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

(60) 
(61) 
(62) 
(63) 


Let us analyze each contributions coming from each terms of Eq. dSU. The first integral Ii given by Eq. (P|l is divergent 
in i?i. In general we can show that 


1 


/ d^r'dz '-—^ 

Iri [42)^+ 4] 


finite n < 2 
00 n > 2 


(64) 


Using the above result we can see that the integrals I3, I4 and the first integral of In are divergent. These integrals 
contain bulk divergences which must be removed introducing bulk counterterms. Next let us analyze the contribution 
coming from the integral I 2 in the region i? 2 - Using previous arguments and considering the external legs we get 


d^rdzd^r'dz' 


IR 2 


K4)" + 4]^ 


< 


d^wdzdz' 


/-e'VO ^0 


(65) 


Therefore this term gives a finite contribution to the Eq. isni). The contribution from the integral Jg to the integral 
must be studied in region i? 2 - In this case we have to consider the external legs, and the property: let us take a 
function f{x,y) positive which does not have singularities except for {x,y) = (0,0), I = Jq/q dxdy f{x,y) then, 
I < fo dx dy fix, y) + 4 dy dy /(x, y), we get 


d^rdzd^r'dz' 


IR 2 


[4-2)^ + 4][(4)^ + 4]^ 


<2 d\ 


dz / du 


z{z + u) 


0 (u^+w^)(u^+z^+w^)^' 


( 66 ) 


From above arguments we have that the contribution from the integral Iq is smaller than 


Sij d^r 

S 4 S 5 
^ S2S3 




ds 


9 


" / + 2 


0 (52+z2)2 

(^ 5 )^ '■ 


+ 


2 S 3 f d^r f dz f du f ds 
J Jo Jo Jo 


zus 


0 (m^ + S^)(m 2 + + z^)^ 


(^ 2)253 


d^r / dw —r < 

Jo 


fSiSo 

\S 2 S 3 ' 


‘{ 82 )^ S 3 


d'^r. 


(67) 
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We conclude that the integral /g is finite. Also integrating the contribution coming from the term /§ on i ?2 we get 

1 


In. [i42r+rh?[{‘^-42r+rU 


( 68 ) 


finite 


Using the fact that the integral d^rdzd^r'dz' j 2 i® finite, we have that this integral also is convergent in 

i? 2 . The contribution from the term Iiq on R 2 is given by 


d^rdzd^r'dz' 


IR2 


[(U2)^ + rl2] M 2 ? + M [(2 - ^2)" + M 


finite 


(69) 


Since the integral U d^rdzd?r'dz' ,, 1 — 5 -r is finite, then the integral defined by Ea. llHnil is convergent 

^ ^ [{z^2y+rf^][{z+^y+rf^] ’ ^ ^ ^ ^ 

in i? 2 . The contribution coming from the terms In contain only a bulk divergence. Otherwise, the contributions 

coming from the terms I 12 is finite. We conclude that we need only bulk counterterms to render the integral defined 

by Eq. (15011 finite. The same analysis can be done for the four-point Schwinger function in the two-loop approximation. 

We obtained that only bulk divergences appear in the full four-point function. 


V. CONCLUSIONS 

In this paper we are interested to show how to implement the renormalization procedure in systems where the 
translational invariance is broken by the presence of macroscopic structures. Eor the sake of simplicity we are 
studying the self-interacting massless scalar field theory in a four-dimensional Euclidean space. We impose that one 
coordinate is defined in a compact domain, introducing two parallel mirrors where we are assuming Dirichlet-Dirichlet 
boundaries conditions. Note that although there are some similarities with the finite temperature field theory using 
the Matsubara formalism, in thermal systems appears only bulk divergences, as for example in the case of system 
where we assume periodic boundary conditions. In non-translational invariant systems, in general to render the theory 
finite it is necessary to introduce surface countertems. 

In this work we generalize some results obtained by Fosco and Svaiter H and also by Caicedo and Svaiter M 
We identify the divergences of the Schwinger functions in the massless self-interacting scalar field theory up to the 
two-loop approximation. First, analyzing the full two and four-point Schwinger functions at the one-loop level, we 
shown that the bulk counterterms are sufficient to render the theory finite. Second, at the two-loop level, we have to 
introduce surface counterterms in the bare lagrangian in order to make finite the full two and also four-point Schwinger 
functions. The most interesting behavior appears in the “double scoop” diagram given by Eq. lim . The amputated 
diagram is given by Eq. lITTI) and we are interested in the last term of Eq. (ESI . This unrenormalized quantity contains 
only bulk divercences. Nevertheless, after the introduction of a bulk counterterm to render the contribution finite 
between the plates, surface divergences appear. Thus this surface divergences must be renormalized. Therefore this 
term contains an overlapping between bulk and surface couterterms. This procedure can be generalized to the n-loop 
level. The inclusion of the counterterm in the lagrangian up two-loop level with the full renormalized action and the 
general algorithm to identify the surface and bulk conterterms in the n-loop level will be left to a future work. 
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APPENDIX A 


In this appendix we will derive an useful representation for the free two-point Schwinger function. Starting from 
Ea. lHTIl . we have that Gg^^(ri — r 2 , Zi, Z 2 ) is given by 


- r2,Zi,Z2) = 


(27r)° 




oip.(ri-r2) 


[p"^ + ' 


(A.l) 


Using the variables u = and v = defined respectively in the region u G [—1,1] and v g [0,2], and also 

making use of a trigonometric identity and performing the sum that appears in Eq. we obtain 


- r 2 ,zi,Z 2 ) = - 




gip.(ri-f2) 

cosh ^A(l — M )(p^ -1- TO^) 2 

cosh ^L(l — z;)(p^-1-m^) 2 ^ 

(p2 + ?n2)i 

sinh -|- 77 j,2)^^ 

sinh ^ A (p 2 -(- 77 i 2 ) ^ ^ 


(A.2) 


(‘ 2 '\ 

Taking m = 0, d = 4, and integrating the angular part, it is possible to show that Gq {ri — r 2 , zi, Z 2 ) can be written 


as 


- f 2 ,Zl,Z 2 ) = 


2(27r)2r'L2 


[ dx( e“’'' - 

cosh ((1 — M ) i) cosh ((1 — n) x) 

/o ^ ^ 

sinh X sinh x 


(A.3) 


where the variable r' is defined by: r' = Making use of the following integral representation of the product 


between the Gamma function and the Riemann zeta function |21| 

r(z)^^ /3 


dx 


'0 


eP^ - 1 


^ —1-1)) 

p 


(A.4) 


where Re{z) > 1, Re{^) > —1 and the Riemann zeta function ^(z, q) is defined by 


OO ^ 


(A.5) 




/o') 

then, it is possible to write Gq Rri — ^ 2 , zi, Z 2 ) as 

1 


- r 2 ,zi,Z 2 ) = 




E 


1 




Finally, using the following identity: 


E 


TT sinh( 27 rr) 


(fc-|)2 + (^)2 


(A.6) 


k— 


(k-zY + r^ 2r sinh^(7rr)-I-sin^(7rz) ’ 


(A.7) 


we obtain the expression for the two-point Schwinger function that we need to proceed in our analysis. Using the 
above equation in Eq. EH) we get 


^(2).^ ^ X sinh(7rr') 


a„(!!a.)sin(2!a) 


[sinh2(2^) -h sin^(2^)] [sinh2(^) -h sin2(2^)] 


(A.8) 
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